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Abstract 

We calculate the electromagnetic form factors of the nucleon within the light-cone sum rule ap- 
proach. In comparison to previous work ID we suggest to use a pure isospin-1/2 interpolating 
field for the nucleon, since the Chernyak-Zhitnitsky current leads to numerically large, unphysical, 
isospin violating contributions. The leading-order sum rules are derived for the form factors and 
the results are confronted with the experimental data. Our approach tends to favor the nucleon 
distribution amplitudes that are not far from the asymptotic shape. 
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1. The elastic scattering of electrons off nucleons at momentum transfer —Q'^ is described by 
the famous Rosenbluth formula IS 



1 + r 2 



(1) 



where (^^(Q^) and Gm{Q'^) are the electric and magnetic Sachs form factors, r = (5^/(4m^), 
m is the nucleon mass and 9 is the scattering angle in the laboratory frame, {da/d^l)-^^^^ is the 
Mott-cross section, which describes the scattering of a pointlike particle. The normalization of the 
form factors at = is given by the nucleon charges and magnetic moments (in units of the 
nuclear magneton, fiN = e/2mp): 

Proton : ^^^(0) = 1, Ga/(0) = /ip = 2.792847337(29) jSj 
Neutron : G's(O) = 0, Gm(0) = /i„ = -1.91304272(45) jSj (2) 

In the Breit frame (^^(Q^) and Gm{Q'^) can be interpreted as the Fourier transforms of the charge 
distribution and magnetization density in the nucleon, respectively. The matrix element of the 
electromagnetic current 0^™(a:) = euu{x)'y ^u{x) + e4d{x)'y ^d{x)) taken between two nucleon 
states is conventionally written in terms of the Dirac and Pauli form factors Fi{Q'^) and F2{Q'^), 
respectively. 



{P-q\jr{0)\P) = N{P-q) 



N{P), (3) 



where P^ is the four-momentum in the initial nucleon state, m is the nucleon mass, P'^ = (P — 
qY = m?, q^ is the (outgoing) photon momentum, Q"^ = —q^, a^^ = |[7^, 7,^] and N{P) is the 
spinor of the nucleon. The electric and magnetic Sachs form factors are related to the Dirac and 
Pauli form factors in the following way 

G'm(Q') = Fi(g') + F2(Q'), Ge{Q^) = F^{Q^) - -^F2(g2). (4) 

It is known that the experimental data for Gm{Q'^) at values of up to 5 GeV^ are very well 
described by the famous dipole formula both for the proton ||4l |5l |6l |71 [HI |9| and for the neutron 
lIT0l[TTl[T2| (following IIT3II we compare our theoretical predictions only with data sets where both 
forward and backward angle data were taken in the same apparatus). 

-G^m{Q') - -Gm{Q') - 2,2 = Gd{Q') ; ^^l - 0.71 GeV^ . (5) 

/in l^i -|- / /igj 

For the electric form factor of the proton the experimental situation currently is unclear. Older 
measurements based on the Rosenbluth separation showed a dipole behavior |l5l |6l |71 [HI |9| of the 
electric Sachs form factor, but in recent measurements at the Jefferson Lab Hall A Collaboration 
using the recoil polarization technique a significant deviation from the dipole was observed |[T?[[T51 
[T^ . This experimental discrepancy has been attracting lots of attention and has not been settled 
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yet (for a review see ifTBlD . The values of the electric form factor of the neutron are very small 

CQiiniiiHi- 

The ultimate goal of the theoretical and experimental analysis of the form factors of the nucleon is 
the determination of the nucleon wave functions. In recent years it has been becoming increasingly 
clear that the strict perturbative approach based on QCD factorization and involving at least two 
hard gluon exchanges is not applicable in the several GeV region and it has to be complemented 
by some non-perturbative techniques. The method of light-cone sum rules (LCSR) lfT9ll suggests 
itself since it incorporates both the perturbative and non-perturbative end-point contributions and 
allows to calculate the form factors as a systematic expansion in terms of nucleon distribution 
functions of increasing twist ll20l 121112211 . Alternative models to determine the form factors of the 
nucleon can be found e.g. in ll23ll . The general concept of the LCSR calculation is familiar from 
numerous applications of this technique to meson decays ll24l and the particular realization for 
baryons was worked out in Ref. ifll. The starting point of the LCSR approach is that one of the 
participating nucleons is substituted by a suitable local current. The choice of the current is a subtle 
issue and is motivated by the necessity to have a strong "nucleon signal" and small sensitivity to 
the contributions of higher resonances and the continuum. In addition, the choice is influenced by 
the particular tasks of the calculation. In particular, in ^ the so-called Chernyak-Zhitnitsky (CZ) 
nucleon current was used since it allows to enhance contributions to the sum rule that are due to 
the leading twist distribution amplitude of interest and suppress higher-twist contributions. The 
main point of this letter is to point out that the use of the CZ current induces large implicit isospin 
violations in the sum rules of order 20% (and more) but this deficiency can be overcome by using a 
modified current which is a pure isospin- 1/2 state. In addition to exact isospin symmetry, using the 
improved current one gets a better stability of the sum rules and a surprisingly good agreement with 
the experimental data using the set of asymptotic distribution amplitudes. We, therefore, argue that 
using the pure isospin current is advantageous and allows to increase the accuracy and reliability 
of the sum rules. Further applications e.g. to axial form factors will be considered in a subsequent 
publication ll25ll . 

2. We start with the electromagnetic coupling of protons and consider the following correlation 
function 

Tr(P,g) =^ J d'xe'''-^{0\T{7^^{0)fJ^ix)}\P) (6) 

which includes an interpolating proton field 77^. The basic principle of sum rules is to calculate 
this correlation function in two ways and finally compare the two results. First one can insert a 
complete set of states between rj^ and j"'^ in Eq. 

TriP, 9) = E ^'(0) ^ - ^) ^2 _ L _ ^1 ^Tio) \p) , (7) 

where A characterizes the state and s stands for the polarization. In lH) the CZ current ll2T]l 

^cz(O) = e''' [u'{0)C A^(0)] 75 ^d'{0) (8) 

was used for 7]^. In this case 

(OIV^IP) = h{Pz)^N{P) (9) 
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(here 2 is a light-cone vector, = 0), and the coupling determines the normalization of the 
leading twist proton distribution amplitude ll2m . Using the definition of the form factors in Eq. Q 
the contribution of the nucleon intermediate state in the correlation function Eq. Q is readily 
derived to be 

^"UP^q) = -^^{P'z){\2F,{Q') (P'z) - F,{Q'){qz) 
m/ — F' I L 



mN,{P), (10) 
m J 



where P' = P — q. In order to get rid of terms ~ Zu that give subdominant contributions on the 
light-cone and to simplify the Lorentz structure we contracted the correlation function with . 
Alternatively, one can calculate the correlation function in Eq. ^ at large Euclidean momenta P'^ 
and = —Q"^ in terms of nucleon distribution amplitudes. To the leading order in the strong 
coupling one gets expressions of the form (cf. Q) 

z-T,{P,q)^i j d^x j ^e^('?+'=)^^(0|e*^V(«ia:)^f(a2x)d2(a3x)|P)a/3^ (11) 

where Cap-y are certain coefficients (involving Lorentz structures) and the real numbers are 
either one or zero. By assumption ~ 1/(P — g)^ and in this limit the remaining three- 
quark operator sandwiched between the proton state and the vacuum can be written in terms of the 
three-quark nucleon distribution amplitudes of different twist t = 3, 4, 5, 6, see Il20ll2ni22t . 



(0| e'^''u'^{a^xy.{a2x)dl{a^x) \P) = ^ F^X^'^F^ (12) 

i 

where F^*^ = V^''\ A*^*), T*^*) are vector, axial-vector and tensor distribution amplitudes and X"^ 
and are Dirac structures which are listed in lE^ . Equating Eq. (fTOt and the QCD calculation at 
a certain intermediate momentum (P — g)^ ~ —1 GeV^ yields a sum rule for the form factors in 
terms of the nucleon distribution amplitudes. The matching procedure involves several technical 
steps that are common for the QCD sum rule approach in general and have the purpose of suppress- 
ing contributions both from higher resonances and the continuum, and of higher-twist operators. 
In particular a Borel transformation is performed, introducing the Borel parameter b instead of 
(P — qY, and the nucleon contribution is defined by introducing a cutoff in the spectral density at 
So ~ (1.5 GeV)^ which is approximately the mass of the Roper resonance. The Borel parameter 
Mb is chosen to be in the range 1.0 — 1.5 GeV, see ||T1|24| for details. 

The nucleon distribution amplitudes that provide the necessary non-perturbative input to the sum 
rules are usually written in terms of the conformal expansion lEHEH. The so-called asymptotic 
distribution amplitudes correspond to taking into account the lowest conformal spin only and com- 
paring the sum rule results with the experimental data one may hope to get an estimate for the 
corrections. In Ref. Q it is shown that large contributions of higher conformal spins are not wel- 
come by the data (the fact that higher terms of the conformal expansion tend to overestimate the 
physical result is already known from the pion form factor [27]), but further work is needed in 
order to make this conclusion quantitative. 
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The question that we address in this work is whether the accuracy of the sum rules can be 
improved by the choice of the nucleon interpolation current. In particular, we look at the isospin 
symmetry. The CZ current (EJ) does not have a definite isospin so that isospin relations between 
different nucleon distribution amplitudes are imposed as the relations between the corresponding 
matrix elements. This current has been chosen for the sum rules in yj because with this choice 
(and in contrast to, e.g. the so-called loffe current ll28ll ) the coefficients C^/j^ in (fTTT) are of order 
one for the contributions of leading twist distribution amplitudes and are suppressed, generically, 
by a power of Af^ for higher twists (to leading order in the strong coupling). The price to pay 
is, however, that in the sum Q there are contributions of both isospin-1/2 and isospin-3/2 states, 
e.g. the A-resonance. It is usually believed that the isospin separation is not important since 
isospin-3/2 resonances are separated from the nucleon by a relatively large mass gap and, therefore, 
sufficiently strongly suppressed by the Borel transformation. One may also speculate that summing 
over states with different isospin in fact makes the spectral density more smooth and thus improves 
the duality approximation for the continuum. Our starting observation is that these arguments can 
be checked by studying the isospin relations for the sum rule predictions. If one determines only 
the electromagnetic form factors of the nucleon, as it was done in lU), the necessity to fulfill isospin 
symmetry is hidden. If, however, one determines in addition to Ff^ (proton in the initial and final 
state) and (neutron in the initial and final state) the form factors F"2' which arise in the vector 
part of the weak-current (j^^'^'^(x) = u{x)^ui^ — l5)d{x)) triggering the /5-decay, one can show 
that the isospin relation 

pnp^ppp_pnn fQ^^ = l2 (13) 

has to hold. Checking whether Eq. (fT3t holds numerically for the sum rule predictions, we can test 
the assumption that the contamination by isospin-3/2 contributions in the sum rules is negligible. 
The corresponding calculations (see |E9) yield the following result: if one uses asymptotic distri- 
bution amplitudes, then the isospin sum rule in Eq. (fTlt is violated by ~ 20 %. If higher conformal 
spin contributions of the distribution amplitudes are taken into account, the isospin violations be- 
come even larger. In other words, the use of the CZ current rjcz for the evaluation of the nucleon 
form factors leads to an unphysical uncertainty of at least 20 % induced by the "pollution" of sum 
rules by the isospin-3/2 contributions. 

The problem can be overcome in a rather simple way by using a modified current which is a 
isospin-1/2 eigenstate. In particular, we suggest to use 

rfj{x) = '^-e'^^ [{u\x)Civ?{x)\ ^4d^{x) - {u\x)C id? {x)\ ^^iu^{x)) (14) 

which is an isospin-1/2 eigenstate and it projects on the leading-twist distribution amplitudes as 
well so that all "good" properties of the CZ current are retained. The factor 2/3 in Eq. ([T?b is 
introduced to fulfill the same normalization condition ©, so that the "hadronic" part of the sum 
rule (fTOb remains intact. On the other hand, using the improved current 77/ for the quark level 
calculation the isospin relations in Eq. (fT3t are recovered exactly. In order to be able to argue that 
the modified current in (fT?b is indeed superiour for the LCSR calculations, we still need to check 
what happens with the sum rule predictions. Since in fH it was found that large corrections to the 
asymptotic distribution amplitudes seem to be in contradiction to the data, in this letter we only 
consider asymptotic distributions as an example. A general case will be studied in (El). The final 
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LCSRs using the improved current 77/ read 
/ 1 

dxi 



2e„ 



3/. 



N 



V? 



2 4 
TTi Tfi 

Pi + Jp^(P^- P3) + ^P4 



{xi) EXPi 



?7i \ / n\ 



Ml 



d X? P4(a;i) 



dxi Q"^ + xf 



g2 + (a;0)2 m? 



EXPo 




jxi - 


^ X3 




/ 1 


2 

m 




Ml 


V- 



'dxi 

Xi 



P2 



m 
Ml 



\ 



m 



P4 (xi) EXPi 
d xi pi{xi) 



I 



dxi Q"^ + xf 



2 n 



Q2 + (a;0)2 ^2 



EXP, 



where for asymptotic distribution amplitudes 



EXPi := exp 



EXP2 := exp 



1-xi Q 
xi 

So — 



m 

2 ^^'m| ' ' 



pi(x) = 60(1 -xf xf. 



N, 



— (1 - x)^ [6x (1 - 4x) Ai + (36 - 370x + 1006x2 - 117x^) /n] 
18 



P3(a:) = - 7^ (1 - ^ [8 (9Ai - 2A2) - 3 (565 - 417x) /, 



P4{X) 



1 - xf x^ [48Ai - 5 (343 - 15x) / 



X, 



180 
1 
2m 



N 



{Q^ + So - m^f + 4m2Q2 _ (q2 ^ _ ^2^ 



(15) 



(16) 



(17) 



The final result depends on the two ratios Xi/ /n and X2/ /n of the non-perturbative parameters 
fN = (5.3 ± 0.5) ■ 10"=^ GeV^, Ai = -(2.7 ± 0.9) ■ 10"^ GeV^ and A2 = (5.1 ± 1.9) • 10"^ GeV^, 
which are discussed e.g. in ll22ll . 

3. The comparison of the sum rule results (IT?!) . (fT^ with the experimental data is shown in 
Figs. 1-5. In all cases the central value of the LCSR prediction is shown by the solid curve while 
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dashed curves show the effect of the variation of the normalization Ai/Zat in the range —5.1 ± 1.7 
which is representative of the possible uncertainty. Varying the Borel parameter Mb in the range 
of 1.2 GeV to 1.6 GeV yielded no sizeable effect; in the plots Mb = GeV is used. 
In Fig. 1 we plotted the magnetic form factor of the proton normalized to the dipole formula. In 
this case the difference compared to using the CZ current appears to be small and our results are 
close to H]. In both calculations the LCSR prediction using asymptotic distribution amplitudes 
tends to overestimate the form factor by about 50%. This disagreement may signal that contribu- 
tions of higher conformal spin have to be included, but in order to make quantitative statements 
one first has to calculate perturbative 0{aj) corrections to the sum rules which is beyond the tasks 
of this letter. The ratio of the electric and the magnetic proton form factors is shown in Fig. 2. Here 
the LCSR prediction is surprisingly close to the experimental values and tends to favor the values 
obtained by the recent experiments at Jefferson Lab lfT4lfT5lfT6l . However, in this case as well, 
without the inclusion of -corrections it is premature to draw definite conclusions. The difference 
to the calculation in fll is quite sizeable for this ratio, up to 50%. In Fig. 3 and Fig. 4 the mag- 
netic and the electric form factors of the neutron are plotted, respectively. The LCSR prediction 
tends to overestimate the magnetic form factor by about 25% while for the electric form factor 
both the experiment and the LCSR give comparable small values. In this cases we again observe a 
noticeable improvement compared to Q. Finally, in Fig. 5 we study the ratio F2/F1 for the proton 
multiplied by Q. We actually plotted QF2/ (KpFi), with the anomalous magnetic moment of pro- 
ton Kp, in order to have the same normalization as the figures in ifT^ . The LCSR calculation shows 
a very weak dependence of this ratio on which agrees with the scaling observed at Jefferson 
Lab lll4lfT5lfT6l . In the LCSR approach such behavior results from an interplay of soft and hard 
contributions with different scale dependence and only holds approximately in a limited range of 
the momentum transfer. 

To summarize, in this work we have presented arguments for the use of the improved nucleon 
current ([T?b in the LCSR calculations. Our current retains all desired properties of the CZ current 
and in addition it is a proper isospin eigenstate so that isospin relations between form factors are 
valid in this case at the level of correlation functions and are reproduced identically in the sum rule 
results for the form factors. Our numerical estimates demonstrate that using the improved current 
one eliminates an implicit uncertainty of the calculations in Q that is due to the isospin symmetry 
violation and also in all cases we obtain a better stability of LCSRs and a better agreement with 
the data using the set of asymptotic three-quark nucleon distribution amplitudes up to twist-6 con- 
structed in II22II . More details and the application to nucleon axial form factors will be considered 
in a forthcoming publication [25 1. It has to be mentioned that the LCSRs to leading-order accuracy 
in the QCD coupling only take into account contributions of "soft" or "end-point" regions that are 
subleading in the true Q"^ ^ 00 limit. The leading contributions appear at the level of perturbative 
corrections to the sum rules and their evaluation presents an important task for further studies. We 
believe that LCSRs with radiative corrections included can provide quantitative information on 
nucleon distribution amplitudes. 

Acknowledgements: We would like to thank V.M. Braun for many enlightening discussions 
and reading the manuscript, N. Mahnke for useful discussions, M.K. Jones for providing the data 
of QF2/F1 which were plotted in fig. 4 of ifT^ and J. Arrington for useful comments on the 
experimental situation. 
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Figure 1: Solid red line: LCSR prediction for the magnetic form factor of the proton normalized 
to the dipole form factor G^/ (/ipGd). Dashed red lines: errors due to the variation of the normal- 
ization Ai//jv. Blue symbols: experimental values: *: SLAC 1994 p9|; A: SLAC 1994 ■: 
SLAC 1970 0*; ♦: Bonn 1971 fTj*; □ : Stanford 1966 |5|* (*: Data actually taken from O)- 
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Figure 2: 5o//J re^? //ne: LCSR prediction for the ratio of the electric and magnetic form factors of 
the proton yUpGE(Q^)/G^(Q^). Dashed red lines: errors due to the variation of the normalization 
Ai/Zat. Red symbols: experimental values obtained via Polarization transfer: ▲: Jefferson LAB 
2002 [[ll'l; ■: Jefferson LAB 2001 O; ♦: Jefferson LAB 2000 ifTi!: Blue symbols: experimental 
values obtained via Rosenbluth separation: ■: SLAC 1994 il; ♦: SLAC 1994 ||9l; A: SLAC 
1970 lEI *; *: Bonn 1971 0*; A: Stanford 1966 O*. (*: Data actually taken from O)- 
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Figure 3: Solid blue line: LCSR prediction for the magnetic form factor of the neutron normalized 
to the dipole form factor G^(Q^) / (AinGD(Q^))- Dashed blue lines: errors due to the variation of 
the normalization Ai //at. Red symbols: experimental values: ■: SLAC 1993 [ilOJ : ▲: Mainz 
2002 03; ♦: Mainz 1998 ED. 
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Figure 4: Solid blue line: LCSR prediction for the electric form factor of the neutron C^iQ"^). 
Dashed blue lines: errors due to the variation of the normalization Xi/ /n- Red symbols: experi- 
mental values: ♦: SLAC 1993 HOI; A: Jefferson Lab 2001 lEI; ■: Mainz 1999 El. 
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Figure 5: Solid red line: LCSR prediction for the ratio QF^Q^) / (KpFf{Q'^)). Dashed red lines: 
errors due to the variation of the normalization \i/ Jn- Red symbols: experimental values obtained 
via Polarization transfer: A, ♦: M. Jones (private communication); Blue symbols: experimental 
values obtained via Rosenbluth separation: ■: SLAC 1994 H; ♦: SLAC 1994 Q. 
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